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I.  INTRODUCTION 


Segmented  penetrators  are  being  investigated  because  of  th<  ir  potential  for  improved 
terminal  ballistic  performance  when  compared  to  monolithic  penetrators.  One  of  the  major 
problems  associated  with  segmented  penetrators  is  how  to  deliver  the  segmented  penetrator 
to  the  target  without  significantly  degrading  the  improved  terminal  ballistic  performance. 
The  Scorpion  projectile  was  developed  to  investigate  one  possible  concept  for  delivering 
segmented  penetrators  to  a  target  under  sponsorship  by  the  Defense  Advanced  Research 
Projects  Agency  (DARPA).  The  concept  involves  a  parent  vehicle  with  segments  trailing  in 
the  wake  of  the  parent  projectile.  The  work  documented  in  this  report  is  concerned  with 
the  aerodynamic  design  of  the  Scorpion  parent  projectile.  This  work  was  also  sponsored 
by  DARPA  and  v;as  carried  out  in  consultation  with  Science  and  Technology  Associates 
(STA),  Inc.,  who  had  overall  responsibility  for  the  projectile  design.  These  results  were 
previously  documented  in  two  Ballistic  Research  Laboratory  (BRL)  interim  memorandum 
reports  (BRL-IMR-939,  March  1990,  and  BRL-1MR-957,  February  1991)  which  were  for¬ 
warded  to  STA,  Inc.  and  DARPA.  In  1992,  the  Scorpion  projectile  was  fired  through  an 
aerodynamics  range  at  the  former  BRL.’  This  final  report  formally  documents  the  work 
discussed  in  the  two  previously  mentioned  interim  reports.  In  addition,  some  experimental 
data  obtained  from  the  aerodynamics  range  firings  are  also  included  to  allow  benchmarking 
of  the  original  predictions. 

Aerodynamic  predictions  for  preliminary  design  configurations,  as  well  as  for  the  cur¬ 
rent  design,  are  documented.  The  predictions  of  the  aerodynamics  of  these  vehicles  were  per¬ 
formed  using  a  sophisticated  computational  technique  which  allowed  the  three-dimensional 
viscous  flow  field  about  the  body  to  be  determined.  Aerodynamic  coefficients  were  deter¬ 
mined  from  the  computed  flow  field.  These  aerodynamic  coefficients  allow  assessment  of  the 
projectile  stability,  pitch-damping  characteristics,  roll  characteristics,  and  in-flight  motion 
(trajectory). 

In  this  report,  the  computational  technique  used  to  determine  the  aerodynamics  of  the 
various  parent  projectile  designs  is  presented  in  the  following  section.  Pitch-plane  aero¬ 
dynamic  predictions  of  two  initial  designs  are  then  discussed.  In  a  subsequent  section, 
pitch-plane  aerodynamics  predictions  for  several  more  recent  design  concepts  (including  the 
current  design)  are  presented  and  compared  with  the  predictions  for  the  initial  designs. 
Comparisons  of  the  predictions  with  the  subsequently  obtained  range  data  are  also  shown. 
Finally,  predictions  of  the  roll  characteristics  of  the  current  design  are  presented. 


II.  COMPUTATIONAL  APPROACH 


Computation  of  the  viscous  flow  field  about  the  flared  projectile  configurations  was 
accomplished  by  solving  the  thin-layer  Navier-Stokes  equations  using  the  parabolized  Navier- 
Stokes  (PNS)  technique  of  Schiff  and  Steger.1  Using  the  parabolized  Navier-Stokes  technique, 

•The  U.S.  Army  Ballistic  Research  Laboratory  was  deactivated  on  30  September  1992  and 
subsequently  became  a  part  of  the  U.S.  Army  Research  Laboratory  (ARL)  on  1  October 
1992. 
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computational  results  are  obtained  by  marching  through  the  grid  from  the  projectile  nose  to 
th  oase.  This  technique  is  applicable  in  the  supersonic  flow  regime  and  requires  that  the  flow 
fi<  contain  no  regions  of  flow  separation  in  the  axial  direction.  Because  the  computational 
approach  requires  only  a  single  sweep  through  the  computational  grid,  it  is  very  efficient 
compared  with  time-marching  approaches  which  require  many  sweeps  through  the  grid. 

The  technique  has  been  applied  successfully  to  a  number  of  projectile  configurations, 
including  axisymrr.:tric  shell,2-3  flared  projectiles,4  and  finned  projectiles.3  Over  the  past 
several  years,  the  computational  technique  has  been  extended  to  allow  predictions  of  the  roll 
characteristics  of  non-axisymmetric  projectiles.5  6  More  recent  extentions  to  the  computa¬ 
tional  approach  permit  predictions  of  the  pitch-damping  aerodynamic  coefficients.7 

As  is  standard  practice,  prediction  of  the  static  aerodynamic  coefficients,  such  as  pitch¬ 
ing  moment  and  normal  force  coefficients,  are  performed  by  computing  the  flow  field  about 
the  projectile  at  a  fixed  angle  of  attack.  The  pitch-damping  coefficients  are  determined  from 
the  side  force  and  moment  acting  on  the  projectile  undergoing  steady  coning  motion.  Steady 
coning  motion  is  defined  as  the  motion  performed  by  a  missile  flying  at  a  constant  angle  with 
respect  to  the  free-stream  velocity  vector  and  undergoing  a  rotation  at  a  constant  angular 
velocity  about  a  line  parallel  to  the  freestream  velocity  vector  and  coincident  with  the  pro¬ 
jectile  center  of  gravity.  This  is  shown  schematically  in  Figure  .  Coning  motion  is,  in  fact, 
a  specific  combination  of  two  orthogonal  planar  pitching  nocior  >,  plus  a  spinning  motion. 
From  linear  flight  mechanics  theory,  it  can  be  shown  that  th  y  itch-damping  coefficients  are 
related  to  the  side  force  and  moment  due  to  steady  coning  motion. 

The  predictions  of  the  roll  characteristics  are  performed  at  CF  angle  of  attack.  The 
roll-producing  moment  can  be  determined  by  computing  the  net  aerodynamic  roll  moment 
acting  on  the  projectile  at  zero  roll  rate.  The  roll-damping  moment,  which  is  by  definition 
the  variation  of  roll  moment  with  roll  rate,  is  determined  by  computing  the  roll  moment  at 
a  number  of  roll  rates.  It  is  interesting  to  note  that  rolling  motion  as  implemented  here  is 
simply  coning  motion  at  0°  angle  of  attack. 

The  flow  field  predictions  of  the  projectile  in  steady  rolling  or  steady  coning  motion  have 
been  performed  in  a  novel  rotating  coordinate  frame  that  rotates  at  the  roll  rate  or  coning 
rate  of  the  projectile.  The  fluid  flow  relative  to  the  rotating  coordinate  frame  does  not  vary 
with  time,  allowing  the  steady  (non-time  varying)  Navier-Stokes  equations  to  be  applied. 
To  implement  the  rotating  coordinate  frame,  t  he  governing  equations  have  been  modified 
to  include  the  effect  of  centrifugal  and  Coriolis  forces.  The  steady,  thin-layer  Navier-Stokes 
equations  are  shown  below. 


d_E_  0F_  dC  -  J_dS_ 

d£  +  dii  +  dC  +  ~  Rc  DC, 


(1) 


Here,  E ,  F,  and  G  are  the  inviscid  flux  vectors,  S  is  the  viscous  flux  vector,  and  //  is  the 
source  term  containing  the  Coriolis  and  centiifugal  force  terms  which  result  from  the  rotating 
coordinate  frame.  Each  of  these  matrices  are  functions  of  the  dependent  variables  represented 
by  the  vector  q{p,pu,  pv,  pw,e),  where  p  and  e  are  the  density  and  the  total  energy  per  unit 
volume,  and  u,  v,  and  w,  are  the  velocity  components  in  x,  j/,  and  2  directions.  The  flux 
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The  pressure,  p,  can  be  related  to  the  dependent  variables  by  applying  the  ideal  gas 


law. 
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The  turbulent  viscosity,  /q,  which  appears  in  the  viscous  matrices,  was  computed  using  the 
Baldwin-Lomax  turbulence  model.8 

The  thin-layer  equations  are  solved  using  the  parabolized  Navier-Stokes  technique  of 
Schiff  and  Steger.1  Following  the  approach  of  Schiff  and  Steger,  the  governing  equations, 
which  have  been  modified  here  to  include  the  Coriolis  and  centrifugal  force  terms,  arc  solved 
using  a  conservative,  approximately  factored,  implicit  finite-difference  numerical  algorithm 
as  formulated  by  Beam  and  Warming.9 

Following  the  approach  of  Schiff  and  Sieger,  the  equations  are  first  linearized  and  placed 
in  delta  form  where  the  equations  are  solved  for  the  difference  in  the  dependent  variables 
rather  than  the  variable  itself.  This  set  of  equations  is  then  factorized  using  the  approach  of 
Beam  and  Warming.  The  following  set  of  equations  is  obtained. 

[A3,  +  (1  -  a)A£  (6„B}  +  DJ)]  Aq‘  =  RHS  (11) 

[*  +  (1  -  cx)A£  -  ±-  (M'/J))]  =  -W  U2) 

-(1  -  a)A${^  [i;:+,(  E/jy  +  ^\F/jy  +  ni+\G/jy } 

+6,  \c"(E/jy  +  C\r/JY  +  ci+l{G/jy] 

+H’-~6<S>)  (13) 

The  form  of  the  equations,  as  well  as  the  notation,  is  similar  to  that  used  by  Schiff 
and  Steger.  Here,  A,  B ,  C,  and  M  are  the  Jacobian  matrices  of  the  flux  vectors  E ,  F,  G, 
and  S.  Further  details  on  the  definitions  of  these  matrices  can  be  found  in  Reference  1. 
The  important  difference  here  is  the  addition  of  the  matrices  D  and  H  due  to  the  rotating 
coordinate  system.  Although  the  Jacobian  matrix,  D,  can  be  included  in  either  the  cir¬ 
cumferential  inversion  or  in  the  normal  inversion,  including  this  term  in  the  circumferential 
inversion  simplifies  slightly  the  implementation  of  the  shock-fitting  boundary  conditions. 

The  computations  presented  here  were  performed  using  a  shock-fitting  procedure  re¬ 
ported  by  Raj  and  t^naussee.10  This  procedure  solves  the  five  Rankine-Hugoniot  jump 
conditions,  two  geometric  shock-propagation  conditions,  and  one  compatibility  equation  to 
determine  the  values  of  the  five  dependent  variables  immediately  behind  the  shock,  as  well 
as  the  position  of  the  shock.  By  including  the  implicit  part  of  the  source  term  due  to  the 
rotating  coordinate  frame  in  the  circumferential  inversion,  the  shock-fitting  procedure  of  Rai 
and  Chaussee  can  be  used  without  modification,  as  long  as  the  correct  free-stream  conditions 
are  specified  as  shown  below  in  nondimensional  form. 

P  =  1 

pu  =  Mucosa  -f  J/Hc  sin  a 
pv  —  Slc{z  cos  a  —  (x  —  xcg)  sin  q  ) 
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The  computational  results  presented  here  were  obtained  using  a  grid  that  consisted  of 
60  points  between  the  body  and  the  shock.  In  the  circumferential  direction,  gridding  was 
performed  over  the  entire  body  (360°).  For  the  axisymmetric  bodies  72  circumferential  grid 
points  were  used.  On  the  finned  or  straked  portions  of  the  projectiles,  300  circumferential 
grid  points  were  used.  In  the  marching  (axial)  direction,  80  marching  planes  were  required 
for  each  caliber  of  body  length. 

The  computations  were  performed  on  a  Cray  X-MP  supercomputer.  For  the  flared 
projectiles,  solutions  for  a  single  set  of  conditions  (Mach  number,  roll  rate,  angle  of  attack, 
etc.)  required  30  to  40  minutes  of  CPU  time.  Solutions  for  the  finned  or  straked  projectiles 
required  about  90  minutes  of  CPU  due  to  the  increased  grid  resolution  required. 


III.  RESULTS  FOR  THE  INITIAL  DESIGNS 


In  this  section,  the  aerodynamic  performance  of  the  initial  designs  is  discussed.  Schemat¬ 
ics  of  the  first  two  design  configurations  are  shown  in  Figures  2  and  3.  Computations  have 
been  performed  to  determine  the  static  and  dynamic  pitch-plane  aerodynamic  coefficients  for 
both  of  the  proposed  flared  projectile  designs.  Using  the  predicted  aerodynamic  coefficients, 
trajectory  simulations  have  been  performed  and  the  pitching  and  swerving  motion  of  the 
projectile  has  been  determined.  The  computations  were  performed  over  a  range  of  Mach 
numbers  ( M ^  —  3.5  to  4.4)  for  free-fiight  (sea-level)  atmospheric  conditions. 

Predictions  of  the  normal  force  and  pitching  moment  coefficient  slopes  were  obtained 
across  the  Mach  number  range  of  interest  for  both  projectiles.  These  predictions  are  shown 
in  Figures  4  and  5.  As  was  mentioned  previously,  projectile  configuration  #  2  has  a  slightly 
longer  cylindrical  body  than  does  configuration  #  1.  This  additional  cylindrical  portion  of 
the  body  does  not  contribute  to  the  lift  of  the  body.  Thus,  the  normal  force  coefficient  slope 
of  both  projectiles  is  nearly  identical.  The  pitching  moment  coefficient  slope  of  projectile  #  2 
is  slightly  higher  than  that  of  projectile  #  1  because  of  the  favorable  location  of  the  center  of 
gravity  relative  to  the  center  of  pressure.  The  predictions  indicate  that  both  configurations 
will  be  statically  stable. 

Figure  6  shows  the  variation  of  the  static  margin  with  Mach  number  for  both  projectiles. 
The  static  margin,  defined  as  the  distance  between  the  center  of  pressure  and  the  center  of 
gravity,  is  normalized  by  the  cylinder  diameter.  The  predictions  of  static  margin  are  in  good 
agreement  with  previous  predictions.11 

The  computed  variations  of  pitch-damping  moment  and  force  coefficients  are  shown  in 
figures  7  and  8.  The  predicted  pitch-damping  force  and  moment  coefficients  are  determined 
from  the  side  fo  ce  and  moment  acting  on  the  projectile  undergoing  steady  coning  motion. 
This  allows  a  steady  flow  computation  to  be  us-'d  to  predict  an  aerodynamic  coefficient  which 
is  usually  associated  with  unsteady  or  time-dependent  motion.  Projectile  #  2  has  a  higher 
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pitch-damping  moment  coefficient  becan.se  of  its  larger  length-to-diameter  ratio.  Despite  the 
larger  pitch-damping  moment  coefficient,  both  projectiles  have  similar  damping  rates.  This 
occurs  because  the  larger  pitch-damping  of  projectile  #  2  is  offset  by  its  increased  moment  of 
inertia.  The  pitch-damping  force  coefficient,  as  shown  in  Figure  8,  does  not  have  a  significant 
effect  on  the  motion  of  tffi.  projectile.  However,  this  coefficient  is  required  to  determine  V 
variation  of  the  pitch-damping  moment  coefficient  with  changes  in  center  of  gravity  loca,  .1. 
These  relations  have  been  determined  by  Murphy.12 

Figure  9  shows  the  development  of  the  pitch-damping  moment  coefficient  over  projectile 
configuration  #  1  at  Mach  4.4.  The  figure  shows  that  most  of  the  pitch-damping  moment  is 
produced  by  the  flared  portion  of  the  body. 

Using  the  aerodynamic  coefficient  predictions,  the  motion  of  the  projectile  undergoing 
planar  pitching  motion  was  determined  for  both  projectiles  and  is  shown  in  Figures  10 
and  11.  The  amplitude  of  the  pitching  motion  has  been  scaled  by  the  initial  yaw  level,  o0. 
The  initial  yaw  level,  as  well  as  the  phase  angle  of  the  yawing  motion,  is  not  known.  These 
parameters  are  functions  of  the  gun  system,  individual  gun  tube,  sabot  package,  and  other 
factors.  Because  of  this,  these  plots  should  not  be  used  to  determine  the  range  position  at 
which  zero  yaw  level  occurs.  Experience  has  shown  that  initial  yaw  levels  as  high  as  3°  might 
be  expected.13 

Both  Figures  10  and  11  show  that  the  pitching  motion  damps  to  about  40%  of  the 
initial  yaw  level  at  about  200  m  from  the  gun.  The  damping  of  both  projectiles  is  similar, 
even  though  piojcctile  #  2  has  a  pitch-damping  moment  coefficient  that  is  20%  gieater  than 
projectile  #  1  This  is  because  the  damping  rate.  A,  is  also  a  function  of  the  mass,  m,  and 
moment  of  inertia,  /.  of  the  body,  as  shown  in  equation  15.  The  increased  pitch-damping 
moment  coefficient  of  projectile  #  2  is  offset  by  its  increased  moment  of  inertia.  The  net 
effect  is  a  difference  in  damping  rate  of  less  than  4%. 

A  =  £|^|Ci,„-CI)--j-(Cm.+C„.)|  (15) 

The  pitching  rate  of  the  projectiles  is  similar,  though  at  300  m  the  projectiles  are  out 
of  phase  by  half  a  cycle. 

The  flue*  mting  part  of  the  swerving  motion  due  to  planar  pitching  motion  is  shown  in 
Figures  12  a7  d  13.  These  figures  show  the  motion  of  the  projectile  center  of  gravity  about 
the  mean  flight  path  of  the  projectile.  The  displacement  of  the  projectile  center  of  gravity 
about  the  flight  path  due  to  this  motion  is  less  than  a  tenth  of  a  caliber  for  1°  of  initial  yaw. 
Of  course,  the  amount  of  displacement  is  proportional  to  the  initial  yaw  level. 

The  motion  of  the  center  of  the  projectile  base  is  also  shown  on  these  figures.  The 
motion  of  the  base  can  be  decomposed  into  two  components;  the  motion  of  the  center  of 
gravity,  and  the  motion  of  the  base  relative  to  the  center  of  gravity.  In  planar  pitching 
motion,  the  center  of  gravity  responds  half  a  cycle  out  of  phase  to  the  pitching  motion.  The 
pitching  motion  also  produces  a  deflection  of  the  base  relative  to  the  center  of  gravity.  This 
motion  is  in  the  same  direction  as  the  centei  of  gravity  displacement.  The  sum  of  these 
motions  produces  deflections  of  about  0.15  calibers  per  degree  of  initial  yaw  at  a  position 
200  in  from  the  gun. 
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The  swerving  motion  of  the  projectile  center  of  gravity  and  the  motion  of  the  base 
center  is  similar  for  both  projectiles. 


IV.  PITCH-PLANE  AERODYNAMICS  OF  IMPROVED 

DESIGNS 


Following  the  analysis  of  the  aerodynamic  performance  of  the  two  initial  designs,  three 
additional  design  concepts  were  proposed  to  improve  the  pitch-damping  performance  of  the 
projectile.  Each  of  these  three  designs  had  a  larger  length-to-diameter  ratio  compared  with 
the  two  initial  configurations.  Additionally,  the  cone  angle  of  the  nose  weis  increased  slightly 
over  initial  configurations.  This  was  done  to  move  the  center  of  gravity  and  the  center 
of  pressure  due  to  nose  lift  further  forward.  Both  of  these  changes  should  improve  the 
pitch-damping  performance.  The  forebodies  on  each  of  the  three  design  configurations  were 
identical.  The  tail  sections  of  the  three  iesigns  consisted  respectively  of  a  15°  flare,  a  15° 
flare  with  six  strakes,  and  a  fin  hub  with  six  fins.  The  strake  and  fin  sweep  angles  (angle 
between  the  projectile  axis  and  fin  leading  edge)  was  initially  specified  as  30°.  However,  in 
this  study,  the  sweep  angle  was  allowed  to  vary  between  15°  and  30°  so  that  the  effect  of 
fin  sweep  could  be  considered.  The  three  designs  are  shown  schematically  in  Figures  14-16. 
Aerodynamic  predictions  for  these  design  configurations  were  performed  at  Mach  4.4  for 
free-flight  atmospheric  conditions. 

The  predictions  of  the  normal  force  and  pitching  moment  coefficient  slopes  and  the 
pitch-damping  force  and  moment  coefficients  for  the  three  designs  are  shown  in  Figures 
17  20.  These  results  are  displayed  as  a  function  of  fin  sweep  or  strake  sweep  angle.  The 
effect  of  changes  in  center  of  gravity  due  to  variations  in  the  fin  or  ctrake  sweep  angle  have 
not  been  considered  in  the  aerodynamic  predictions.  These  changes  are  probably  small, 
given  that  the  center  of  gravity  location  of  the  flared  projectile  and  the  flared  projectile  with 
30°  strakes  are  similar. 

The  aerodynamic  coefficients  for  the  flared  projectile  with  strakes  increase  as  the  sweep 
angle  of  the  strakes  increases.  The  increases  in  these  coefficients  is  greatest  between  17.5°  and 
25°  sweep  angles.  There  is  little  change  in  the  aerodynamics  coefficients  between  25°  and  30° 
sweep  angle.  The  aerodynamic  coefficients  for  the  flared  projectile  with  strakes  approaches 
the  values  for  the  flared  projectile  as  the  sweep  angle  is  decreased.  This  is  expected  since 
the  flared  projectile  with  strakes  reduces  to  a  flared  projectile  when  the  strake  angle  is  equal 
to  the  flare  angle  (15°). 

The  finned  projectile  shows  a  small  variation  in  the  aerodynamic  coefficients  with  in¬ 
creasing  sweep  angle.  The  normal  force  and  pitching  moment  slopes  foi  the  finned  projectile 
are  slightly  greater  than  those  for  the  flared  projectile.  The  pitch-damping  coefficients,  how¬ 
ever,  are  significantly  less  than  those  for  either  the  flared  projectile  or  the  flared  projectile 
with  strakes 

The  predicted  pitch-damping  coefficients  have  been  used  to  determine  the  yaw-damping 
as  a  function  of  range  for  each  of  the  three  designs  (Figure  21 ).  The  yaw-damping  predictions 
are  compared  with  the  predictions  for  initial  design  projectile#  1.  The  flared  projectile  with 


strakes  has  the  best  yaw-damping  performance.  The  yaw  level  at  200  m  from  the  gun  is  half 
that  of  the  initial  design.  The  flared  projectile  shows  some  improvement  in  yaw-damping 
over  the  initial  design,  while  the  finned  projectile  shows  only  a  marginal  improvement.  It 
should  be  noted,  however,  that  the  finned  projectile  has  considerably  less  drag  than  any  of 
the  flared  projectile  designs.  This  makes  fins  more  suitable  for  tactical  configurations  which 
must  maintain  their  velocity  over  longer  ranges  than  are  required  for  the  Scorpion  pi  jectile. 

In  the  current  application,  one  of  the  primary  benefits  of  the  flare  is  the  increase  in 
the  size  of  the  wake  region.  It  is  believed  that  the  increased  size  of  the  wake  region  may 
decrease  the  aerodynamic  interaction  of  the  trailing  segments.  Because  the  finned  configura¬ 
tion  doesn’t  provide  any  benefit  in  pitch-damping  compared  with  the  flared  projectiles  and 
because  the  finned  projectiles  will  have  a  smaller  and  possibly  more  complex  wake  region, 
they  are  probably  not  suitable  for  the  current  application. 

Based  on  these  aerodynamic  predictions  discussed  here,  STA,  Inc.  has  chosen  to  proceed 
with  the  flared  projectile  with  25°  strakes.  The  25°  strakes  yield  almost  identical  performance 
to  the  30°  straked  projectile  and  is  less  likely  to  encounter  interference  problems  due  to  sabot 
discard. 

The  variation  in  the  normal  force,  pitching  moment,  and  pitch-damping  coefficients 
with  Mach  number  for  the  flared  projectile  with  25°  strakes  is  shown  in  Figures  22-25. 
The  coefficients  show  moderate  variations  between  Mach  3  and  Mach  5.  The  pitching  and 
swerving  motion  of  this  design  are  shown  in  Figures  26  and  27.  These  figures  show  that  the 
pitching  motion  and  the  swerving  motion  at  the  center  of  gravity  and  projectile  base  have 
been  reduced  compared  with  the  initial  design  configurations.  Due  to  an  increase  in  the 
pitching  moment  coefficient,  the  projectile  also  has  somewhat  higher  pitching  rate  compared 
with  the  initial  designs. 

Predictions  of  the  forebody  drag  have  been  made  with  the  parabolized  Navier-Stokes 
(PNS)  code.  This  .ode  is  not  suitable  for  computing  the  flow  in  the  base  region  of  the 
projectile;  hence,  the  base  drag  could  not  be  determined  using  this  approach.  Estimates  of 
the  base  drag  were  obtained  using  the  Naval  Surface  Weapons  Center  (NSWC)  aerodynamic 
design  code.14  The  base  drag  component  was  added  to  the  PNS  forebody  drag  predictions  to 
obtain  predictions  of  the  total  zero  yaw  drag  of  the  projectile.  These  predictions  are  shown 
in  Figure  28.  The  increase  in  drag  for  this  configuration  over  the  initial  configurations  is  due 
primarily  to  the  steeper  flare.  Increasing  the  flaie  angle  from  12°  to  15°  added  about  30% 
to  the  forebody  drag.  The  addition  of  the  strakes  to  the  flare  added  about  5%  additional 
forcbody  drag. 


V.  ROLL  CHARACTERISTICS  OF  CURRENT  DESIGN 


Statically  stable  projectiles,  such  as  the  configurations  examined  here,  arc  designed  to 
fly  at  a  small  roll  rate.  This  is  done  to  minimize  the  effect  of  asymmetries  which  may  occur 
during  manufacturing.  Computations  were  performed  to  determine  the  roll  characteristics 
of  the  flared  projectile  with  25°  swept  stakes.  The  strakes  were  canted  to  generate  a  roll- 
producing  moment.  Figure  29  is  a  schematic  showing  the  nature  of  the  canted  strakes.  The 
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computations  were  performed  over  a  range  of  supersonic  Mach  numbers  and  0°  angle  of 
attack  for  free-fiight  atmospheric  conditions. 

Computations  were  first  performed  to  predict  the  effect  of  cant  angle  on  the  roll- 
producing  moment.  Figure  30  shows  the  variation  in  the  roll-producing  moment  with  cant 
angle  at  Mach  4.4.  The  variation  of  the  roll-producing  moment  shows  a  slightly  nonlinear 
variation  with  cant  angle.  The  largest  nonlinear  variations  were  observed  at  the  higher  Mach 
numbers  examined  here.  The  curve  is  well-fit  using  the  nonlinear  form  shown  below. 

Clo  =  Clt6  +  C,t7\6\6  (16) 

The  coefficients  Ci6  and  Cit7  were  determined  using  a  least-squares  fitting  procedure.  These 
coefficients  are  shown  in  Table  1.  Attempts  were  made  to  fit  the  data  with  a  cubic  term 
instead  of  the  quadratic  term,  however,  the  quadratic  term  provided  a  better  fit. 


Table  1.  Roll-Producing  Moment  Coefficients  as  a  Function  of  Mach  Number 


MACH 

NUMBER 

CU 

cn(J 

3.0 

5. IS 

7.32 

3.5 

5.96 

1.45 

4.0 

6.13 

-1.91 

4.4 

5.98 

24.18 

5.0 

G.C1 

27.45 

The  roll-damping  moment,  Cip,  was  determined  by  running  a  series  of  computations  at 
various  spin  rates.  The  variation  in  the  net  roll  moment  with  spin  rate  at  Mach  4.4  is  shown 
in  Figure  31 .  Computational  results  are  shown  for  three  different  cant  angles,  0.25°,  0.5°,  and 
0.625°.  The  variation  of  the  net  roll  moment  with  spin  rate  (by  definition  the  roll-damping 
moment),  is  seen  to  be  linear  and  independent  of  cant  angle.  For  Mach  4.4,  the  predicted 
roll-damping  moment  coefficient  is  8.95.  The  variation  in  the  roll-damping  coefficient  with 
Mach  number  is  shown  in  Figure  32. 

Figure  31  also  shows  that  the  equilibrium  spin  rate  (spin  rate  at  w-hich  the  net  roll 
moment  is  zero)  increases  with  the  cant  angle.  The  variation  in  the  equilibrium  spin  rate 
with  Mach  number  for  the  0.5°  canted  strake  is  shown  in  Figure  33.  Also  shown  here  is 
the  variation  of  the  nondimensional  pitching  frequency  with  Mach  number.  The  equilibrium 
spin  rate  for  the  0.5°  canted  strakes  is  below  the  pitching  frequency  across  the  Mach  number 
range. 

F  igure  34  shows  the  predicted  roll  history  of  the  projectile  for  the  0.5°  canted  strakes, 
launched  at  Mach  4.4  and  zero  initial  spin  rate.  The  effect  of  the  variation  in  the  projectile 
velocity  due  to  drag  has  been  included  in  the  prediction.  Also  shown  is  the  yawing  frequency 
of  the  projectile.  The  yawing  frequency  varies  slightly  with  range  due  to  the  change  in  the 
projectile  velocity  due  to  drag.  The  spin  history  of  the  projectile  is  well  below  the  pitching 
frequency  over  the  first  200  m  of  flight. 


VI.  COMPARISONS  BETWEEN  PREDICTIONS  AND 

FREE  FLIGHT  DATA 


In  1992,  a  number  of  the  25°  straked  flare  projectiles  were  fired  through  BRL’s  Transonic 
Range.  The  primary  purpose  of  the  tests  was  to  examine  the  segmentation  process  as  it 
occurred  in  flight.  However,  the  projectile  did  fly  over  100  m  before  the  segmentation  process 
was  initiated.  Standard  aerodynamic  range  measurements  were  obtained  throughout  the 
flight  and,  from  this  data,  it  was  possible  to  extract  the  aerodynamic  forces  and  moments. 
To  avoid  disturbances  from  the  initiation  of  the  segmentation  process,  only  measurements 
from  the  first  three  groups  of  shadowgraph  stations  were  used  for  extracting  the  aerodynamic 
data  on  the  parent  projectile.  Additionally,  many  of  the  rounds  had  very  low  yaw  and 
little  useful  data  on  the  aerodynamics  of  the  parent  projectile  was  obtained.  It  should  be 
emphasized  that  the  test  was  primarily  to  observe  the  segmentation  process  and,  to  this  end, 
the  projectile  was  designed  to  fly  with  a  minimum  amount  of  yaw. 

Most  of  the  rounds  yielded  good  data  on  the  drag  coefficient  of  the  round.  Comparison 
of  the  predicted  and  measured  drag  coefficients  are  shown  in  Figure  35.  Both  results  show  a 
similar  trend  with  Mach  number  with  the  prediction  about  5%  lower  than  the  measurement. 
This  is  a  reasonable  comparison  given  that  predicted  base  drag  was  based  on  estimates  from 
a  fast  design  code. 

Comparisons  of  pitching  moment,  pitch-damping  moment,  and  normal  force  coefficients 
with  the  range  data  are  made  in  Figures  36-38.  Since  the  pitching  moment  and  pitch 
damping  moment  are  dependent  on  the  center  of  gravity  location,  the  results  also  include 
predictions  of  these  coefficients  at  the  center  of  gravity  location  as  measured  before  firing. 
The  center  of  gravity  location  as  tested  was  0.27  calibers  forward  of  the  original  design 
center  of  gravity  location.  The  new  aerodynamic  predictions  were  obtained  using  the  center 
of  gravity  translation  relations  and  the  original  base  line  aerodynamic  coefficients.  The 
predictions  are  generally  within  the  scatter  of  the  range  data.  The  accuracy  of  the  range 
data  for  these  coefficients  was  diminished  somewhat  due  to  the  low  yaw  experienced  by  the 
projectile  and  due  to  the  shortened  trajectory  over  which  aerodynamic  data  on  the  parent 
projectile  could  be  obtained.  Again,  it  is  emphasized  that  extracting  aerodynamic  data  on 
the  parent  projectile  was  a  secondary  consideration  for  the  test. 


VII.  CONCLUSION 


Aerodynamic  coefficient  predictions  have  been  obtained  for  several  proposed  designs 
for  the  Scorpion  parent  projectile.  These  aerodynamic  coefficient  predictions  were  used  to 
assess  the  stability  and  damping  performance  of  the  designs.  The  predictions  show  that  the 
flared  projectile  with  25&  strakes  yielded  better  damping  performance  compared  to  the  initial 
designs. 

Predictions  of  the  roll  characteristics  of  the  flared  projectile  with  25°  strakes  were  also 
made.  The  results  show  that  canting  the  strakes  at  0.5°  with  respect  to  the  projectile  axis 
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will  produce  a  roll  rate  which  is  below  the  pitching  frequency  of  the  projectile  during  the 
first  250  m  of  flight. 
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ALL  DIMENSIONS  IN  CALIBERS  (ONE  CALIBER  -  40  6  mm) 
Figure  2.  Schematic  of  projectile  configuration  #  1. 


ALL  DIMENSIONS  IN  CALIBERS  (ONE  CALIBER  -  40  0  mm) 
Figure  3.  Schematic  of  projectile  configuration  #  2. 
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MACH  NUMBER 


Figure  4.  Variation  of  normal  force  coefficient  slope  with  Mach  number. 
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Figure  5.  Variation  of  pitching  moment  coefficient  slope  with  Mach  number. 
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Figure  0.  Variation  of  static  margin  with  Mach  number. 


MACH  NUMBER 

Figure  7.  Variation  of  pitch-damping  moment  coefficient  with  Mach  number. 
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Figure  8.  Variation  of  pitch-damping  force  coefficient  with  Mach  number. 
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Figure  9.  Development  of  pitch-damping  moment  coefficient  over  projectile  #  1,  Mach 
4  4 
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Figure  10.  Pitching  motion  as  a  function  of  range  -  Projectile  #  1. 
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Figure  11.  Pitching  motion  as  a  function  of  range  -  Projectile  #  2. 


y 


16 


DISPLACEMENT  (CALIBERS)  «  DISPLACEMENT  (CALIBERS) 


0.15 


PER  DEGREE  OF  YAW 


-0.10 


-0.15 


1  \ 

■  '  /  \  A 

/'N  \  /  *  /  i 


0.00 


0.05  'V' 


1  /  \  /  \ 


U  H  w  V 

\  !  i  /  \  / 


I  \ 


'  I  1  / 

V- -  MOTION  OF  CG 

—  MOTION  OF  BASE  CENTER 

.  | 

100  200 

RANGE  -  METERS 


ire  12.  Swerving  motion  as  a  function  of  range  -  Projectile  #  1. 
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Figure  13.  Swerving  motion  as  a  function  of  range  -  Projectile  #  2. 
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ALL  DIMENSIONS  IN  CALIBERS  (ONE  CALIBER  =  40  6  mm) 
Figure  14.  Schematic  of  flared  projectile  configuration. 


ALL  DIMENSIONS  IN  CALIBERS  (ONE  CALIBER  «  40  6  nr) 
Figure  15.  Schematic  of  flared  projectile  with  strakes. 


ALL  DIMENSIONS  IN  CALIBERS  (ONE  CALIBER  «  40  6  mm) 
Figure  16.  Schematic  of  finned  projectile  configuration. 
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Figure  17.  Variation  of  normal  force  coefficient  slope  with  sweep  angle. 
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Figure  18.  Variation  of  pitching  moment  coefficient,  slope  with  sweep  angle. 
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Figure  19.  Variation  of  pitch-damping  moment  coefficient  with  sweep  angle. 
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Figure  20.  Variation  of  pitch  damping  force  coefficient  with  sweep  angle. 
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RANGE  -  METERS 


Figure  21.  Yawing-damping  as  a  function  of  range  for  the  improved  designs  compared  with 
one  of  the  initial  designs. 
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Figure  22.  Variation  of  normal  force  coefficient  slope  with  Mach  number,  flared  projectile 
with  25“  strakes. 
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Figure  23.  Variation  of  pitching  moment  coefficient  slope  with  Mach  number,  flared  pro¬ 
jectile  with  25°  strakes. 
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Figure  24.  Variation  of  pitch-damping  moment  coefficient  with  Mach  number,  flared  pro¬ 
jectile  with  25°  strakes. 
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Figure  25.  Variation  of  pitch  damping  force  coefficient  with  Mach  number,  flared  projectile 
with  25°  strakes. 
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Figure  26.  Pitching  motion  as  a  function  of  range,  flared  projectile  with  25°  strakes. 
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RANGE  -  METERS 

Figure  27.  Swerving  motion  as  a  function  of  range,  flared  projectile  with  25°  strakes. 
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Variation  of  drag  coefficient  with  ivlach  number,  flared  projectile  with  25° 
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Figure  29.  Schematic  showing  canted  strakes. 


Figure  30.  Variation  of  roll-producing  moment  coefficient  with  cant  angle. 
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Figure  31.  Variation  of  roll  moment  with  spin  rate  for  various  cant  angles. 
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Figure  32.  Variation  of  roll-damping  coefficient  with  Mach  number. 
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igure  33.  Variation  of  equilibrium  spin  rate  with  Mach  number. 
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Figure  34.  Spin  history  for  0.5°  cant  angle. 
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Figure  36.  Comparison  of  measured  and  predicted  pitching  moment  coefficient,  flared 
projectile  with  25°  strakes. 
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Figure  37.  Comparison  of  measured  and  predicted  pitch-damping  moment  coefficient,  flared 
projectile  with  25°  strakes. 
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Figure  38.  Comparison  of  measured  and  predicted  norma!  force  coefficient,  flared  projectile 
with  25°  strakes. 
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LIST  OF  SYMBOLS 


Cd 

Cl, 

Cto 

C:, 


Cl 


c., 

cLa 

Cma 

cm„  +  cm, 

CV,  +  C.vo 

D 

e 

E,F}G 

H 

1 

J 

m 

A/oo 

P 

pD/V 


Poo 

Be 

S 

$rej 

U,V,W 

u,v,w 

x,y,z 


freestream  speed  of  sound 

drag  coefficient 

net  roll  moment  coefficient 

roll-producing  moment  coefficient 

variation  of  roll-prcducing  moment  coefficient 

with  cant  angle,  6 

variation  of  roll-producing  moment  coefficient 
with  square  of  cant  angle 
roll-damping  coefficient 

slope  of  the  lift  coefficient  with  angle  of  attack 

slope  of  the  pitching  moment  coefficient  with  angle  of  attack 

pitch-damping  moment  coefficient 

slope  of  the  normal  force  coefficient  with  angle  of  attack 
pitch-damping  force  coefficient 
projectile  diameter 

total  energy  per  unit  volume,  nondimensionalized  by  px 

flux  vectors  in  transformed  coordinates 

source  term  resulting  from  rotating  coordinate  frame 

projectile  transverse  moment  of  inertia 

Jacobian 

projectile  mass 

freestream  Mach  number 

spin  rate,  as  used  in  roll  equations  and  roll  coefficients 

riondimensional  spin  rate,  nondimensionalized  by  body 

diameter,  D,  and  free  stream  velocity,  V' 

pressure,  as  used  in  thin-layer  Navier-Stokes  equations, 

nondimensionalized  by  poo 

freestream  static  pressure 

nondb  .  isionalized  by 

Reynolds  number,  VooPooD/ Pco 

viscous  flux  vector  in  transformed  coordinates 

reference  cross-sectional  area  of  projectile,  7rZ)2/4 

velocity  components  in  x,  y,  and  z  directions, 

nondimensionalized  by  a 

Contravariant  velocities  of  the  transformed  Navier-Stokes  equations 
Cartesian  coordinates  with  respect  to  the  body,  nondimensionalized  by  D 


Mote:  Moment  coefficients  are  scaled  as  follows;  A//^/>0ua£oA/£>ZXS're/ 
Force  coefficients  are  scaled  as  follows;  F/ ^Po0a'l0M%0Srcf 


Greek  Symbols 
a  angle  of  attack 

7  ratio  oi  specific  heats,  as  used  in  Navier-Stokes  equations 
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6 

cant  angle  of  strakes 

A 

damping  rate  of  projectile 

P 

laminar  viscosity 

Pt 

turbulent  viscosity 

transformed  coordinates  in  Navier-Stokes  equations 

p 

density 

poo 

freestream  density 

Slc 

coning  rate  of  projectile,  nondimensionalized  by  £>/ 
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